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agreement with the predictions of the self-consistent theory for the 2D Anderson localization. 
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The metal-insulator Anderson transition plays a cen¬ 
tral role in the study of quantum disordered systems. 
Using a tight-binding description of an electron in a lat¬ 
tice, Anderson [1] postulated in 1958 that the dominant 
effect of impurities in a crystal is to randomize the di¬ 
agonal term of the Hamiltonian, and showed that this 
may lead to a localization of the wavefunction, in sharp 
contrast with the Bloch-wave solution for a perfect crys¬ 
tal. In a weakly disordered (3D) crystal, the eigenstates 
are delocalized, leading to a diffusive (metallic) trans¬ 
port, while strong disorder produces an insulator with 
localized eigenstates. From its original solid-state physics 
scope [1-4] this approach has been applied to a large class 
of systems in which waves propagate in disorder. This 
includes quantum-chaotic systems [5, 6] and electromag¬ 
netic radiation [7-9]. Important theoretical progress was 
obtained in ref. [10], which postulated that Anderson lo¬ 
calization can be described by a one-parameter scaling 
law, leading to the prediction that, for d < 2, the dy¬ 
namics is generically localized, even if the disorder is very 
weak. For d > 2, it predicted the existence of the An¬ 
derson transition between a diffusive dynamics at weak 
disorder and a localized dynamics at strong disorder. 

There is no fully quantitative theory of Anderson lo¬ 
calization, and analytic results are scarce. Supersym¬ 
metry techniques [11] allow derivation of expansions in 
powers of d — 2 of the various quantities of interest, but 
reaching even d = 3 is difficult. A useful, simplified the¬ 
oretical approach is the so-called self-consistent theory 
of localization. In few words, it can be thought as a 
mean field theory where large fluctuations are neglected, 
but where weak localization corrections to transport, due 
to interference between time reversed multiply scattered 
paths, are included self-consistently. For spinless time- 
invariant systems, belonging to the orthogonal symmetry 
class [11], this approach correctly predicts the existence 
of the metal-insulator Anderson transition for d > 2, al¬ 
though it fails to predict the correct critical exponent. 


For d = 1, it quantitatively predicts the localization 
length in a weak disorder. Other approaches lead to ap¬ 
proximate values for the critical exponent not far from 
the numerical prediction [12]. 

Dimension d = 2 - the lower critical dimension - is 
very special, the localization properties depending on the 
symmetry class. In the orthogonal symmetry class, the 
dynamics is always localized, but the localization length 
is predicted to scale exponentially with the inverse of the 
disorder strength, i.e. ^ oc ^e^x.p{'Kk^/2) [13] where k is 
the wavevector and i the mean-free path for propagation 
in the disordered medium. As discussed in the Supple¬ 
mental Material [14], such an exponential dependence is 
a signature of the fact that d = 2 is the lower critical 
dimension for Anderson localization. The 2D case has 
been previously studied experimentally in optical and ul- 
tracold atom systems [8, 15], but no quantitative indi¬ 
cation of the exponential scaling has been demonstrated 
yet. In the present Letter, we use the well-known corre¬ 
spondence between the d-dimension Anderson model and 
the d-frequency quasiperiodic kicked rotor [6, 16, 17] to 
test experimentally these predictions. 

The quasiperiodic kicked rotor (QPKR) [5, 6, 16-18] 
is a spatially one-dimensional system with an engineered 
time-dependence such that its dynamics is similar to the 
dynamics of a time-independent multidimensional sys¬ 
tem. The QPKR can be simply realized experimentally 
by exposing laser-cooled atoms (Cesium in the present 
work) to a delta-pulsed (kicked) laser standing wave of 
wavenumber kL and time period Ti. The amplitude of 
the kicks is quasiperiodically time-modulated with a fre¬ 
quency co’ 2 . The dynamics is effectively one-dimensional 
along the axis of the laser beam, as transverse directions 
are uncoupled. The corresponding Hamiltonian is: 

2 ^-1 
H = -—\- K cos x[l e cos {uj 2 t)] S{t — n) , (1) 

where x is the particle position, p its momentum, K the 
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kick intensity and 5 the amplitude of the modulation. We 
have chosen conveniently scaled variables such that dis¬ 
tances along the x axis are measured in units of 
the particle’s mass is unity and time is measured in units 
of the pulse period Ti. In the quantum case, a crucial pa¬ 
rameter is ^ = Ahkj^Ti/M^ the reduced Planck constant. 

For 5 = 0, one obtains the periodic kicked rotor, which 
can be mapped onto a one-dimensional Anderson-like 
model [6], and displays “dynamical” localization [5, 19], 
that is, Anderson localization in momentum space in¬ 
stead of configuration space. For non-zero 5, the tem¬ 
poral dynamics of the QPKR is exactly that of a two- 
dimensional periodic kicked system [16, 17, 20], which it¬ 
self can be mapped - provided 27r/Ti, uj 2 and k are incom¬ 
mensurable numbers - onto a two-dimensional Anderson- 
like anisotropic-hopping model, where anisotropy is con¬ 
trolled by 5 and the ratio of hopping to diagonal disorder 
is controlled by K/k [21]. 

The existence of a mapping of the kicked rotor onto an 
Anderson-like model has been used to experimentally ob¬ 
serve ID Anderson localization with atomic matter waves 
as early as 1995 [18]. The two-frequency modulation of 
the QPKR - which can be mapped on a 3D Anderson 
model [17, 20] - was used to experimentally observe 3D 
Anderson localization and the metal-insulator Anderson 
transition [22], measure accurately the critical exponent 
and demonstrate its universality [23]. 

The experimental study of the 2D case is more chal¬ 
lenging than the 3D one, because the observation of the 
exponential behavior of the localization length pioc re¬ 
quires pioc to vary over about one order of magnitude. 
The localization time increasing with pioc, this also re¬ 
quires the ability to preserve coherence over several hun¬ 
dreds kicks. This needed major evolutions of our experi¬ 
mental setup [24]. 

Experimentally, an atomic sample consisting of few 
million atoms is prepared in a thermal state (3.2 /iK) 
whose momentum distribution is much narrower than the 
expected localization length. The atomic cloud is then 
“kicked” by a far-detuned (A ^ 13 GHz) pulsed standing 
wave (SW). Pulse duration is typically r = 300 ns, while 
the typical pulse period Ti = 27.778 /rs corresponds to an 
effective Planck constant k = 2.89. According to Eq. (1), 
an adjustable amplitude modulation with C(;2/27r = a/ 5 
is superimposed to the kick sequence. In our previous 
experiments, to minimize coupling with gravity, the SW 
was horizontal. However, for 1000 kicks the atoms fall 
down by 3.8 mm, compared to the 1.5 mm SW waist, lim¬ 
iting the maximum number of kicks to 200. In order to 
overcome this limit we used in the present experiment a 
vertical SW, and the atoms fall freely between kicks. The 
new SW is formed by two beams that can be indepen¬ 
dently controlled, both in amplitude and phase, through 
a radio-frequency driving two acousto-optic modulators. 
This allows us to accurately cancel gravity effects, by 
imposing a linear frequency chirp to one arm of the SW 


with respect to the other, so that the SW itself “falls” 
with acceleration g. A kicked rotor is thus realized in 
the free-falling reference frame. These technical improve¬ 
ments are discussed in more detail in the Supplemental 
Material [24]. At the end of the sequence, the veloc¬ 
ity distribution is measured by a standard time-of-flight 
technique. 

Eigure 1(a) shows experimental momentum distribu¬ 
tions n(p) recorded after 0 to 1000 kicks for K = 
5.34,^ = 2.89,5 = 0.36. If the dynamics were classical, 
the momentum distribution would keep its initial Gaus¬ 
sian shape and the average kinetic energy increase lin¬ 
early with time, = Ekin(^ = 0) + Dt, where D is 
the classical diffusion constant in momentum space. In 
contrast, the experimental result displays a distribution 
which diffusively broadens at short times, but tends to 
freeze, i.e. to localize at long times. This clear-cut proof 
of localization is confirmed by the shape of the momen¬ 
tum distribution, shown in Eig. 1(b) after 200 kicks. It 
very clearly displays an exponential shape [25] (a straight 
line in the logarithmic plot) exp (— |p|/pioc) /2pioc charac¬ 
teristic of localization with a localization length pioc [26]. 
Eigure 1(c) shows the momentum distributions after 1000 
kicks for K = 5.34, ^ = 2.89 and increasing values of 5. 
It demonstrates that the localization length varies very 
rapidly with 5, indicating the evolution from a ID local¬ 
ization at 5 = 0 to a truly 2D localization with a much 
longer localization length at 5 = 0.6. In order to pre¬ 
vent trivial localization on KAM tori [27], we always used 
> 4, ensuring that the classical dynamics is ergodic. 

Instead of measuring the full momentum distribution, 
it is sufficient to measure the population Hq (t) of the zero 
velocity class as 


“ 4n»(t) 

is proportional to (p^)(t) (as the total number of atoms 
is constant) [28]. 

Eigure 2 displays (at 1000 kicks) vs. 5 for var¬ 
ious values of K and k, showing that the exponential 
dependence in 5 is a general feature, with a rate that 
decreases with k and increases with K. Note the overall 
^kin dynamics of a factor of 60 (corresponding to a 8 fold 
increase in the localization length), a key feature of the 
present experiment. 

The scaling theory of localization [10] predicts that di¬ 
mension d = 2 is the lower critical dimension for the 
Anderson transition. Eor a time-reversal invariant spin¬ 
less system (thus belonging to the orthogonal universality 
class), all states are localized with an exponentially large 
localization length. Eor a usual 2D time-independent sys¬ 
tem, the relevant parameter is the dimensionless conduc¬ 
tance at short scale, equal to the product ki of the wave 
vector by the mean free path, so that the logarithm of 
the localization length is proportional to kl [14]. 
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Figure 1. (Color online) Experimentally recorded momentum 
distributions for the kicked rotor exposed to a quasi-periodic 
driving, Eq. (1). (a) K = 5.34,^ = 2.89,£ = 0.36, 0 to 

1000 kicks (step 200). The momentum distribution diffusively 
broadens at short times and freezes at longer times, proving 
the existence of 2D Anderson localization. Time increases 
from top to bottom curves, (b) Momentum distribution at 
200 kicks in log scale, showing the exponential shape charac¬ 
teristic of localization. The circles are experimental points, 
the blue dashed line is a Gaussian fit and the black solid 
line an exponential fit for \p\ > 3(2/i/cl). (c) Localized mo¬ 
mentum distributions after 1000 kicks, as a function of the 
anisotropy parameter £, for K = 5.34, k = 2.89 as in (a) 
and (b). The modulation amplitude £ increases from top to 
bottom curves. The rapid increase of the localization length 
shows the evolution from the ID localization at £ = 0 to the 
truly 2D Anderson localization. Note the different horizontal 
scales in the various plots. 
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Eigure 2. (Color online) Kinetic energy Ekin of the quasiperi- 
odic kicked rotor vs. the modulation amplitude £, for various 
values of the kicking strength K and effective Planck constant 
k. The error bars indicate the typical experimental uncer¬ 
tainty. The four curves are straight lines in this logarithmic 
scale, with a slope that decreases with k and increases with 
K. 


The scaling theory cannot be directly transposed to the 
case of the kicked rotor for two reasons: i) There is no 
wavevector playing the role of k. Instead, one must con¬ 
sider the diffusion constant (in momentum space), which 
is, for a periodic kicked rotor, approximately equal to 
iF^/4. ii) The diffusion process for the 2D quasiperiodic 
kicked rotor is not isotropic. As shown in [29] and dis¬ 
cussed in the Supplemental Material [14], the quasiperi¬ 
odic kicked rotor can be mapped on a 2D Anderson-like 
model, whose dynamics at short time is indeed diffusive, 
but anisotropic. Along the “physical” direction (which co¬ 
incides with the atom momentum component along the 
standing wave), the diffusion constant is - for small 5 
- almost equal to the one of the periodic kicked rotor. 
Dll ~ along the other (virtual) direction, the dif¬ 

fusion constant is D 22 ~ K‘^£^ so that it vanishes in 
the limit e ^ 0, where one must recover the usual ID 
periodic kicked rotor. 

Altogether, the relevant parameter is the geometric av¬ 
erage of the diffusion constant along the two directions 
\/DiiD 22 oc eK^. The scaling theory predicts that the 
logarithm of the localization length should be propor¬ 
tional to ^/DiiD^/k^. A similar prediction was made 
in [16] using a slightly different method. 

The self-consistent theory of localization is an attempt 
towards more quantitative predictions, based on the same 
ideas as the scaling theory. It has been successfully used 
to predict properties of the Anderson transition [30, 31], 
and was transposed to the periodic kicked rotor in [19, 32] 
and to the quasiperiodic kicked rotor with two addi¬ 
tional driving frequencies in [33]. It consists in comput¬ 
ing perturbatively the weak localization correction to the 
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Figure 3. (Color online) Increase in the kinetic energy at 
t = 1000 (oc Pi^oc) of fho quasiperiodic kicked rotor with re¬ 
spect to the purely one-dimensional situation £ = 0 vs. the 
scaling parameter The cloud of experimental points 

- collected at various values of K, s and ^ - is distributed 
around an average linear dependence in this semi-logarithmic 
plot, which shows the exponential dependence of the localiza¬ 
tion length, characteristic of 2D Anderson localization. The 
red dashed line is the prediction of Eq. (3). The spread is due 
in part to experimental imperfections (at large 5(iF/^)^, the 
localization time is not much shorter than the duration of the 
experiment) and in part to fundamental reasons: The linear 
dependence on sK^ in the argument of the exponential, 
Eq. (3), is valid only at small £, and the formula assumes that 
the classical diffusion constant is proportional to ^ while 
the actual diffusion constant has oscillatory corrections. The 
black curves are numerical simulations corresponding to the 
two “extreme” values of Kj^ — 1.3,^ = 3.46 (higher curve) 
and Kj^ — 2.5,^ = 2.89 (lower curve); they display the same 
spreading phenomenon. 


(anisotropic) diffusion constant and to extrapolate to the 
strong localization regime. It however depends on the 
cut-offs used. For our quasiperiodic kicked rotor [14] it 
confirms the prediction of the scaling theory, namely: 


^ioc = ^exp 


(olsK^ 


( 3 ) 


where a is a number of the order unity, which may how¬ 
ever depend smoothly on the parameters. In the limit 
5 ^ 0, it is a = tt/v/32- 

In Fig. 3, we display the results of 275 measurements, 
corresponding to 12 values of the ratio Kj^ G [1.3, 2.5], 
with K G [4.33,7.26] and ^ = {2.89,3.2,3.46}, and to 
5 values from 0 to 0.6 (step 0.06). Dividing F^kin(^) by 
E\cin{s = 0) makes it possible to probe the exponential 
term in Eq. (3). The exponential dependence (straight 
line in logarithmic scale) is visible for 5 < 1, materialized 
by the red dashed line, corresponding to the prediction 


a = of the self-consistent theory. Despite the 

spreading of the experimental results around the average 
trend, the overall agreement is rather good. This proves 
the exponential dependence of the localization length in 
2D, and thus that d = 2 is the lower critical dimension 
for the metal-insulator Anderson transition. Some devi¬ 
ations are nevertheless visible. They arise from different 
phenomena: First, for large 5, the localization time can 
be only slightly shorter than the duration of the experi¬ 
ment (1000 kicks), meaning that the measured momen¬ 
tum distribution is not the asymptotic one for infinite 
time and underestimates the eventual saturation of F^kin 
at long time. This explains why the experimental points 
at large 5 tend to lie below the theoretical prediction. 
This is confirmed by numerical calculations in the ex¬ 
perimental conditions for the largest value of Kj^ = 2.5 
(longest localization time), see the solid lower curve in 
Fig. 3. A second, more fundamental, phenomenon is that 
Eq. (3) assumes that the classical diffusion constant is 
simply iF^/4, which is valid only for iF ^ 1, whereas 
oscillatory corrections at moderate K are known to exist 
for the ID kicked rotor [34] and to persist even for the 3D 
QPKR [29]. This dependence is thus not eliminated by 
the normalization to E’kin(^ = 0). This explains a signifi¬ 
cant part of the spreading of the data. Finally, Eq. (3) is 
expected to be valid in the e ^ 0 limit, see Supplemen¬ 
tal Material [14]. At larger 5 values, higher order terms 
must come into play and are responsible for significant 
deviations. This is visible in Fig. 3, where both experi¬ 
mental (points) and numerical (solid lines) data are well 
predicted at small eK^ ^ but are more widely spread 
as increases. A thorough analysis of all these de¬ 

viations is beyond the scope of this Letter. 

To summarize, we presented the first experimental 
evidence of two-dimensional Anderson localization with 
atomic matter waves. We studied the variation of the lo¬ 
calization length with the system parameters and showed 
that it displays an exponential dependence characteris¬ 
tic of time-reversal spinless systems. To the best of our 
knowledge, such an experimental evidence has not been 
observed previously. It demonstrates experimentally that 
d = 2 is the lower critical dimension of the Anderson 
transition. The observed localization length varies as pre¬ 
dicted by the scaling and the self-consistent theories of 
localization. 
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Appendix: Supplemental Material 

The following Supplemental Material discusses: 

I. The main improvements in our experimental setup 
allowing the 2D Anderson localization to be ob¬ 
served. 

II. The mapping of the quasiperiodic kicked rotor onto 
a 2D Anderson model, indicating the essentials of 
this mapping and its peculiarities relevant for the 
present work. 

III. The classical diffusion in the 2-frequency quasiperi¬ 
odic kicked rotor. 

IV. A schematic derivation of the localization length, 
indicating the essential steps of the calculation in 
the case of the quasiperiodic kicked rotor with 2 
incommensurate frequencies and the origin of the 
exponential dependence on the anisotropy. 

1. Setup improvements 

In our previous experiments, in order to minimize cou¬ 
pling with gravity, the standing wave (SW) was horizon¬ 
tal, see Eig. 4a. This geometry has several drawbacks. 
Eor 1000 kicks, the atoms would fall of 3.8 mm, to be 
compared to the 1.5 mm of the SW waist. Moreover, 
the SW was built by retro-reflecting the incoming beam, 
which leads to a ^I m optical path difference between 
the two beams overlapping on the atomic cloud region, 
a source of phase noise detrimental to dynamical local¬ 
ization. These restrictions limited the kick number to 
less than 200 in our previous experiments. In order to 
overcome this limit, we built a new SW system with sev¬ 
eral improvements, as shown in Eig. 4b. The SW is now 
vertical, and, between kicks, the atoms fall freely. The 
intensity inhomogeneity in the transverse direction was 
also reduced. The SW is now formed by two indepen¬ 
dent beams, which has many advantages, as each arm 
can be independently controlled, both in amplitude and 
phase through the radio-frequency wave that drives the 
acousto-optic modulators. This allows us to accurately 
cancel gravity effects, by imposing a linear frequency 
chirp to one of the arms with respect to the other, so that 
the SW itself “falls” with acceleration g. A kicked rotor 
is thus realized in the free-falling reference frame. Ei- 
nally, the SW phase noise induced by the laser linewidth 



Figure 4. Evolution of the experimental setup, a) Horizontal 
standing wave setup used in former experiments. The kicks 
are controlled by a single acousto-optic modulator (AOM), 
and the standing wave is built by retro-reflecting the incom¬ 
ing laser beam, b) Schematic view of the new vertical setup 
used in the present work. A Distributed Feedback Laser Diode 
(DFB) seeds two optical amplifiers. The kick-sequence tem¬ 
poral modulation is produced by an arbitrary wave generator 
(AWG). In addition to the pulsed sequence, a linear chirp is 
added so that the standing wave “falls” with acceleration g, 
simultaneously with the atomic cloud. 

is minimized by accurately balancing the optical paths 
of the arms to better than I cm. This is performed by 
directly minimizing the kinetic energy dispersion at 1000 
kicks with a ID kicked rotor. 

2. Mapping of the quasiperiodic kicked rotor onto 
a 2D Anderson model 

The mapping of the quasiperiodic kicked rotor 
(QPKR) onto a 2D Anderson model has been described 
in [16, 17, 20]. Because it is useful to derive an approxi¬ 
mate expression for the localization length, we summarize 
here the basics of the calculation. 

The starting point is the following Hamiltonian: 

1-L = -^^uj 2P2^K cosxi [1 -h 5 COSX 2 ] ^(5(t-n) (A.I) 

n 

which describes a periodically kicked, two-dimensional 
“rotor”; we use here quotations marks around “rotor” be¬ 
cause of the unusual form of the kinetic energy along the 
direction 2, linear instead of quadratic. It is a matter 
of algebra [20] to show that the temporal evolution of 
an initial state initially localized at X 2 = 02 , but with 
arbitrary wave-function along xi\ 

T(xi, X2, t = 0) = E(xi, t = 0)J(X2 - (j)2) (A.2) 
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leads to the state: 

^{xi,X2,t) = E{xi,t)6{x2 - 02 -^2t), (A.3) 

where the wavefunction S(xi,t) exactly obeys the time- 
dependent Schrodinger equation of the quasiperiodic 
kicked rotor: 

H = + K cos xi[l e cos {uj 2 t + 02)] S{t — n) . 

n 

(A.4) 

Thus, the evolution of the quasiperiodic kicked rotor can 
be thought of as the evolution of a 2D periodic “rotor”, 
with a peculiar initial state perfectly localized in the X 2 
direction and thus completely delocalized in the conju¬ 
gate p 2 direction. 

The Hamiltonian Eq. (A.l) being time-periodic, one 
can use the Floquet theorem and look at its Floquet 
eigenstates. Let |0) denote an eigenstate of the evolution 
operator over one period with eigenvalue exp(— 
and consider the state |x) = {1 iW)~^\(f)) where 

W{xi,X 2 ) = tan [K cos xi{l e cos X 2 )/ 2k]. The Hamil¬ 
tonian being periodic in xi and X 2 (we take periodic 
boundary conditions), |x) can be expanded on the basis 
of plane waves |mi,m 2 ) - eigenstates of the momentum 
operators in both directions with eigenvalues rriik and 
lx) = Y.rm,m2Xmt,m2\m.i,m2). It is again a mat¬ 
ter of algebra [6, 17] to show that the coefficients Xmi,m 2 
obey the following equation: 


• The hopping is anisotropic, governed by the K/k 
coefficient along direction 1 and by the eK/k co¬ 
efficient along direction 2. For 5 = 0, the 2D sys¬ 
tem appears as a series of chains along direction 1, 
which are uncoupled along direction 2, and a ID 
Anderson model is recovered. 

• Floquet quasi-eigenstates with different quasi¬ 
energies E are associated with the same energy 0 
in Eq. (A.5), but with different realizations of the 
disorder, Eq. (A.6), all having the same statistical 
properties. This is why, for given values of param¬ 
eters LC, 5, all quasi-eigenstates have the same lo¬ 
calization length [14]. 

In order to understand localization properties of the 
quasiperiodic kicked rotor, it is thus sufficient to study 
transport and localization on the equivalent anisotropic 
Anderson-like 2D model Eq. (A.l). 


3. Classical diffusion 

It is well known that the classical dynamics of the peri¬ 
odic kicked rotor is described by the Standard Map [27]. 
A similar map can be constructed for the quasiperiodic 
kicked rotor. The classical evolution over one period of 
the Hamiltonian Eq. (A.4) is given by the following map: 


T ^ ^ ^^ri,r2Xmi+ri,m2-|-r2 ^5 (^•^) 

7’i,r27^0 


where 


-mi ,7712 


is given by: 


= tan < - 


. rni 


. E 

+ Ld2'^2 ) — 


(A.6) 


and the Wr^,r 2 two-fold Fourier components of 

the doubly-periodic function W. 

Equation (A.5) can be interpreted as the eigenvalue 
equation for a Anderson-like model on a 2D lattice with 
sites labeled (mi, m 2 ), with hopping described by HAi,r 2 
and on-site energies €mi,m 2 ' There are however four dif¬ 
ferences with respect to a usual Anderson model: 

• The hopping matrix elements Wr^,r 2 limited 

to nearest neighbors. They do however decrease 
fast enough at large (ri,r 2 ) so that this difference 
does not change the qualitative behavior (i.e. lo¬ 
calization). 


Piu+i =Plr^ + A:sinxi^(l -h5C0SX2^) , 

P2^+1 =P2r. E Kecosxi^smx2^ , 

^In + l = ^In +Pln + 1 ^ 

^2n + l = + ^2 , 


where = y{t = n E e), y = Xi,pi {i = 1,2). If K 
is sufficiently large, the classical dynamics is almost fully 
chaotic [29]. For the ID problem (5 = 0 ), this takes place 
for K > 6. For the 2D problem (non vanishingly small 
5), we found an almost fully chaotic dynamics for iX > 4, 
and this is why we did not perform any experiment below 
this value. In the chaotic regime, the kicks in momentum 
have a pseudo-random sign (depending on xi and ^ 2 ), 
making the dynamics in momentum space appear as a 
pseudo-random walk, i.e. a chaotic diffusive process at 
long time. Note however that the diffusion is anisotropic 
because kicks along p 2 are typically smaller than kicks 
along Pi by a factor 5. In the limit of large K it is easy 
to evaluate the diffusion tensor, defined as 


• The on-site energies emi,m 2 random, but 

rather a deterministic pseudo-random sequence. 
Provided tt, k, UJ 2 are incommensurate, the sequence 
has no periodicity and localization properties sim¬ 
ilar to those of a truly random model are expected 
(and observed numerically). 


Don = lim 

t —^-j-oo 


{Pipj) 

t 


(A.7) 


by assuming that positions of consecutive kicks are uni¬ 
form uncorrelated variables (see [29] for the essentially 
identical calculation in 3D). It is diagonal in the (1,2) 
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directions with: 

Dii«(7fV4)(l + eV2), 

D22 - K^e ‘^/8 , 

Di^j ss 0 . (A.8) 

Numerical simulations of the classical dynamics [29] fully 
confirm that the classical dynamics is an anisotropic dif¬ 
fusion; however, at not too large K, correlations between 
successive kicks are responsible for oscillatory corrections 
to the diffusion constant of the order of few ten percents. 


4. Localization length 


The calculation of the localization length turns out to 
a bit tricky in 2D. It has been discussed in the literature 
mainly in the context of electrons in disordered poten¬ 
tials [35, 36]. The calculation for the quasiperiodic rotor 
follows the same lines. We sketch here only the general 
method, comparing the key steps with the case of elec¬ 
trons in disordered potentials. 

The starting point is to take into account the weak lo¬ 
calization effect due to closed loops in the system. In the 
weak scattering regime, where this correction is small, it 
takes the following form for electrons (or massive parti¬ 
cles) in disordered potentials: 


r 1 dq 
D 7Tp J ( 27 r)^ Dq‘^ 


(A.9) 


where D is the classical (Boltzmann) diffusion constant 
and V the quantum modified one. The term after 1 in 
the right-hand side is the weak localization correction. 
It depends on the density of states p and involves a two- 
dimensional integral over the vector q which is a momen¬ 
tum, conjugate of position. The integrand is nothing but 
the classical diffusive kernel (at zero frequency, hence infi¬ 
nite time). Formally, the integral diverges both for small 
and large g, so that appropriate cut-offs must be set. The 
natural short distance (large q) cut-off is the mean free 
path I (at shorter scale, the dynamics is not diffusive). 
The short q cut-off can be taken proportional to l/L, 
where L is the size of the system. One then obtains a 
size-dependent diffusion constant V{L). Taking for the 
density of states in 2D its disorder-free value l/27r (we 
take the mass of the particle and h as unity), one obtains: 


correct for arbitrary large size, as the correction is com¬ 
puted assuming it is small. This ceases to be true when 
the right hand side term in Eq. (A. 10) vanishes. This 
gives an order of magnitude of the size at which the dif¬ 
fusion constant vanishes, that is the localization length. 
One then gets 


C = (A.ll) 


Of course, this is only a very approximate expression. 
The self-consistent theory of localization [35, 36] is an 
attempt to be a bit more quantitative. It predicts essen¬ 
tially the same exponential dependence for the localiza¬ 
tion length. 

For the periodic kicked rotor, the weak localization 
corrections have been computed in [19]. These results 
have been extended to the quasiperiodic rotor in [33, 37]. 
There are essentially two modifications: 

• The term depending on the density of states (which 
is meaningless for a time-periodic system where the 
density of states of the Floquet Hamiltonian is in¬ 
finite) 1/iTp must be replaced by . 

• Because the classical dynamics is an anisotropic 
diffusion, the diffusive kernel is now {Duql -f 
7^22^2When performing the integral over q, it 
is not entirely clear how to choose the cut-offs. The 
simplest choice is to take the large qi and q 2 cut-offs 
scaling like the (anisotropic) mean free paths, i.e. 
respectively proportional to 1/^/D^ and 1 /a/D^. 
This is however arbitrary and questionable. It is 
important to understand that the choice of differ¬ 
ent cut-offs will affect the prefactors, but not the 
key point, namely the logarithmic dependence of 
the integral on the system size. 


With this simple choice of cut-offs, the weak localization 
correction reads: 


■DiiiL) ^ V22{L) ^ ^ 

Dll D22 'K\/DiiD22 V ^ / 


(A.12) 


where I is a yet unspecified short scale cut-off. 

With the same reasoning as the one for electrons, one 
obtains the approximate expression for the localization 
length (in momentum space for the kicked rotor): 


D{L) 

D 




(A.IO) 


Pioc = I exp 


^\/DiiD22^ 


(A.13) 


where k is the wavevector of the particle. In these units, 
the classical diffusion constant is D = ki/2. The logarith¬ 
mic dependence in L is a crucial ingredient, as we will see 
below. It is intimately related to the low-g divergence of 
the integral. 

Equation (A.IO) shows that the diffusion constant de¬ 
creases with the system size. Of course, this cannot be 


In the 5^0 limit, one must recover the localiza¬ 
tion length of the periodic kicked rotor, which fixes 
I = K‘^ IAk. By inserting the values of the diffusion con¬ 
stants, Eqs. (A.8) into Eq. (A.13), we finally obtain: 


K‘^ 

Pioc = exp 


7rsK‘^ 


(A.14) 





Because of the questionable assumptions on the cut¬ 
offs, it might well be that the coefficient 7r/A/32 is not 
exact and requires some correction. This is a difficult 
question left for further investigation. We just emphasize 
that the exponential behavior of the localization length 
with the scaling parameter comes from the pref¬ 

actors in the integral over q and from the logarithmic 
singularity of the integral. It is thus a very robust phe¬ 
nomenon, unlikely to be affected by the shortcomings of 
the self-consistent theory of localization. Ultimately, the 
experimentally observed exponential dependence of the 
localization length proves that the weak localization cor¬ 
rection is logarithmic with the system size. This, in turn, 
implies that dimension d = 2 is the lowest critical dimen¬ 
sion of the Anderson transition. Indeed, a modification 
of the dimension would introduce an extra factor in 
the integral, i.e. would change the logarithm dependence 
in L in Eq. (A. 12) to an algebraic one, incompatible with 
the exponential behavior of the localization length. 
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